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Spectra and magnetic properties of large spins J, placed 
into a crystal electric field (CEF) of an arbitrary symmetry 
point group, are shown to change drastically when J changes 
by 1/2 or 1. At a fixed field symmetry and configuration of its 
N extrema situated at p-fold symmetry axis, physical charac- 
teristics of the spin depend periodically on J with the period 
equal to p. The problem of the spectrum and eigenstates 
of the large spin J is equivalent to analogous problem for a 
scalar charged particle confined to a sphere S 2 and placed 
into magnetic field of the monopole with the charge J. This 
analogy as well as strong difference between close values of J 
stems from the Berry's phase occurring in the problem. For 
energies close to the extrema of the CEF, the problem can be 
formulated as Harper's equation on the sphere. The 2 J + 1- 
dimensional space of states is splitted into smaller multiplets 
of classically degenerated states. These multiplets in turn are 
splitted into submultiplets of states transforming according 
to specific irreducible representations of the symmetry group 
determined by J and p. We classify possible configurations 
and corresponding spectra. Experimental realizations of large 
spins in a symmetric environment are proposed and physical 
effects observable in these systems are analyzed. 

PACS numbers: 03.65.Sq, 03.65.Bz, 75.10.Dg, 02.20.Df 



I. INTRODUCTION 

Conventional wisdom accepts that large spins or or- 
bital momenta J (in units of K) are almost classical. In 
particular, if J ^> 1, their measurable properties do not 
change substantially if J changes by 1/2 or 1. This com- 
mon belief was undermined by Haldane Q who demon- 
strated that the ground-state and spectrum of the low- 
energy states in one-dimensional spin chains are abso- 
lutely different for integer and half-integer spins. 

In this paper we show that similar phenomena can be 
observed on the level of an individual spin placed into 
external electric field. If the field possesses high symme- 
try (cubic or icosahedral) the distinction between spins 
becomes more subtle. For example, in the case of cubic 
symmetry not only integer spins differ from half-integer 
(this difference is intuitively obvious due to the Kramers 
degeneracy), but the remainder at division of the spin 
by four occurs to determine the spectrum and degener- 
acy of the low-lying states. These striking differences 
can be found in experiment either by spectral analysis 
or by magnetic measurements. We will show that spins 
1000, 1001, and 1002 placed into a cubic environment 



have 100% different magnetic susceptibilities at low tem- 
perature. Moreover, we will show that a kind of random- 
ness appears in properties of large spins in some cases 
and variation of large spins by one can change magnetic 
and spectral properties in incontrollable way. 

Certainly, the conventional wisdom we started with is 
presumably correct. It is wrong only in a very small 
range of energy or temperature, the smaller the larger is 
J. Nevertheless, as it already happened with the Haldane 
theory, these deviations from classical behavior may be 
important for the experiment. 

The source of all these peculiarities is the Berry's 
phase. Physically, it is associated with the fact that, 
when the classical rotator moves on its unit sphere, it si- 
multaneously rotates around its axis. The rotation phase 
distinguishes the rotator from a quantum or classical par- 
ticle confined on a sphere. The rotator problem can 
be reduced to the particle problem, but the represent- 
ing particle must have an electric charge of unity and 
must be subjected to the homogeneous magnetic field of 
a monopole with the magnetic charge J placed into the 
center of the sphere. In quantum mechanics J accepts 
integer and half-integer values. 

This paper is composed as follows. In the next section 
we introduce quasi-classical description of large spins. 
The Berry's phase, Berry's connection, and reduction to 
the problem of a charged particle in the monopole field 
arc considered in Section [II. In the fourth section we 



perform the group analysis of the problem. The fifth sec- 
tion contains the derivation of the low-energy spectrum 
and magnetic properties of large spins. We separated the 
case of random levels in Section [v|. Numerical calcula- 
tions for a special po tent ial in a wid e rang e of spin values 
are given in Section VII. In Section |VIII we propose ex- 
perimental realizations of large spins. Our conclusions 
can be found in Section IX. 



Brief reports on a part of this work were published 
earlier P,|l3[]. 



II. QUASI-CLASSICAL DESCRIPTION OF 
LARGE SPINS 

The classical image of a large spin is the classical ro- 
tator, i.e., a vector with a fixed length J. Its position is 
determined by two spherical coordinates 9 and 4>. Some- 
times coordinates J z = J cos 9 and <f> are more convenient 
since they have a simple Poisson brackets: {J 2 ,0} = 1. 
Classical motion is determined by the Hamiltonian: 
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H = /(J)-h-J, 



(1) 



where h is magnetic field (with a precision of a constant 
factor) and /(J) is an arbitrary function of J, invari- 
ant with respect to inversion: J — * —J. The latter 
requirement is equivalent to the time reversal symme- 
try 0. Together with the standard Poisson brackets 
{Ji,Jj} = EijkJk the Hamiltonian (Q) contains full in- 
formation on classical spin dynamics. Periodical trajec- 
tories on the sphere can be quantized according to the 
Bohr quantization rule: 



Jr. 



,E)d<j> 



7bK 



(2) 



E 



where J z ((ft,E) can be found from equation /(J) 
with the substitution: J x = y/ J 2 — J 2 cos (ft, J y = 
\J J 2 — J 2 sin (ft, and 7b is a constant. 

Let us first consider general properties of spin trajec- 
tories in zero magnetic field. The function /(J), being 
continuous on the sphere, has at least two minima and 
two maxima. If the external crystal field has a non-trivial 
symmetry group, the number of equivalent minima is 
larger. For example, it can be equal to 4 for tetrago- 
nal symmetry, 6 for hexagonal symmetry. In the case of 
cubic symmetry it can be 6, 8, or 12 (directed along 4-, 
3-, and 2-fold axes respectively). The number of equiva- 
lent minima for icosahedral symmetry can be 12, 20, and 
30 (directed along 5-, 3-, and 2-fold axes respectively). 
We considered the situations when extrema are located 
in the symmetrical positions. In principal, it is possible 
that they are in more general asymmetric positions. 

Classical trajectories can be separated into two classes: 
"localized" and "delocalized" . If energy is close enough 
to the minimum (maximum) of /(J), the trajectories are 
confined in a vicinity of one of the minima (maxima). 
We call such trajectories localized. In the intermediate 
region of the energy trajectories are "delocalized" , they 
are not confined near any of the extrema. It is obvious 
that delocalized trajectories are highly model-dependent, 
i.e, they depend on a specific form of /(J). Localized 
trajectories are much more universal: they depend only 
on the symmetry and on the positions of the minima. The 
same remark is correct with respect to quantized levels: 
low- lying levels, close to / m i n , or almost maximal values 
of energy, close to / max , have universal features, whereas 
levels in between are rather non-universal. Therefore, 
further we will study only a part of the spectra close to 
/min or / max . Note that the spectrum of the quantum 
problem is discrete and limited by f mm and / max - 

Before we proceed to detailed study of these levels let 
us make an important remark. For any fixed J and any 
given /(J) the quantum problem consists in the diago- 
nalization of (2 J + 1) x (2 J + 1) matrix. Therefore, the 
question arises whether the general theory is necessary. 
The answer is yes. First of all because no reliable infor- 
mation about function /(J) is available. We present here 
general facts, independent on specific form of /(J), but 



only on its symmetry group and specific configuration 
of the extrema. The only requirement for our theory is 
J> 1. 

Thus, classically a localized stationary state is multi- 
ply (iV-fold) degenerate. Quantum fluctuations provide 
a finite radius for each of these states which can be enu- 
merated as |1), |2), . . . , \N). For considered large J all 
these states are oscillatory ones within the precision 1/ J. 
More subtle, but not least essential, quantum effect is the 
tunneling between these states. The tunneling amplitude 
between two states \i) and \j), i ^ j, is exponentially 



small 



ex exp(— CijJ), where cy are constants for a 



given /(J). Therefore, we take into account only tun- 
neling between the nearest-neighbor states, i.e., the ones 



with the smallest 



c, and neglect tunneling between 



more remote states with > c. To estimate the value 
of c, we need to specify the Hamiltonian. For simplic- 
ity we consider the case of the cubic symmetry with the 
Hamiltonian: 



H? = -a{Jt + J 4 + J 4 ) 



(3) 



where a > is a constant. The minimum value of 
is E m i n = —a,J A . There are six minima corresponding 
to the directions of the 4-fold axes: (± J, 0, 0), (0, ± J, 0), 
(0, 0,±J). Let us consider, for example, tunneling be- 
tween minima (J, 0,0), (0, J, 0). By symmetry the tun- 
neling trajectory is the smaller arc of the big circle pass- 
ing through these points (Fig. [l]). Putting H° — £mm 
we find from eqn. (|): 



Jz(4>) = ±iJ\ 



1 — cos A(ft 
7 + cos A(ft ' 



(4) 



The tunneling amplitude is proportional to the exponent; 



,7T/2 

oc exp(i / J z ((ft)d(ft) = exp(— ( Jln3)/2) = e 
Jo 



For a more realistic Hamiltonian 



-0.55J 



njO 



h° - b{4 + 4 + jl + mjljpl), 



(•5) 



(6) 



the exponential factor in the tunneling amplitude is 
exp(— c(u)J), where c(u) is a function of the ratio u = 
bJ 2 /a. The graph of c(u) is shown in Fig. || for values 
of u in the interval —2/3 < u < 1/15 (a > 0), where 
the tunneling path passes along the geodesies. In the 
region 1/15 < u < 3 (a > 0), the six minima are still 
global, however, the tunneling trajectories (there are two 
of them due to the symmetry) deviate from the geodesies 
(see Fig. ^) and the estimation of the exponent becomes 
more complicated. Effects of the multiple path tunnel- 
ing, for the case of the octahedron configurations, will 
be considered in Section A numerical analysis of the 
Hamiltonian (|^) and comparison to the predicti ons o f the 
semiclassical approximation is given in Section VII. 
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Another important feature of the Hamiltonian (g) is 
that, depending on signs of a and b and the parameter 
u, it displays 6, 8, or 12 minima. The phase diagram for 
this, important for applications, Hamiltonian is shown in 
Fig. |j| On the boundaries of different "phases" different 
groups of minima become equal each other. Then, in the 
quantum problem, the degeneracy of the ground state 
increases, for example, from 6 to 14. Therefore one can 
expect some singularities in close vicinity of the bound- 
aries. 

For the case of the icosahedral symmetry the simplest 
Hamiltonian is: 



i^x i^x ~~ Jy ) Jy i^y 



■ mjljpl ~ 3a/5 X 



(7) 



J 2 Z ) + J 2 Z J 2 X (J 2 Z - J 2 ))), 



The minimum value of TC[ is aJ 6 /5 (a < 0). There are 
12 minima corresponding to the vertices of an icosahe- 
dron (directions of the 5-fold axes): J(a, j3, 0), J(0, a, /?), 
J(/3,0,a), where a 2 = (5+\/5)/10 and f3 2 = (5-V5)/10. 
A calculation similar to the one for TL® gives the ex- 
ponential part of the tunneling amplitude exp(— 0.28 J). 
Positive values of a yield the 20-fold configuration with 
the minima along the 3-fold axes. 

Addition of the next non-trivial invariant of the icosa- 
hedron group, a polynomial of the 10-th order over J, 
allows the configuration with 30 minima along the 2-fold 
axes. 

The tunneling partly lifts the classical degeneracy. 
What was the TV-fold degenerate state without tunneling 
is splitted into a multiplet of sub-levels separated by ex- 
ponentially small energy intervals oc exp(— cJ), whereas 
the distances between different multiplets arc propor- 
tional to 1/ J. Each sub-level in the multiplet corresponds 
to a finite-dimensional subspace of states transforming 
according to an irreducible representation of the symme- 
try group. However, as we mentioned already, the real- 
ization of this group and the spectrum for the rotator is 
very different from those for a quantum particle confined 
on a sphere. Anyway the problem is reduced to diag- 
onalization of a square matrix of the rank TV (classical 
degeneracy of the level) with non-zero matrix elements 
between geometrically closest states only. We neglect 
the tunneling between more remote states (non-nearest- 
neighbors) unless otherwise stated. 




FIG. 1. Tunneling trajectories of the spin (single paths) 
The 6-fold configurations of O. 
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FIG. 2. c(u) for the 6-fold configuration of 
the single tunneling path regime. 
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FIG. 3. Phase diagram of Hamiltonian (g) in bj 2 -a plane. 
The dashed lines separate the regions of the single and double 
path tunneling. The numbers on the periphery are the slopes 
of the corresponding lines. 
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III. BERRY'S PHASE, BERRY'S CONNECTION. 

In the framework of quasi-classical spin dynamics the 
spin is treated as a rigid vector fixed by its direction 
n. The closest quantum analog is the so-called coherent 
state |n) which is defined as an eigenstate of operator n- J 
with the maximal eigenvalue J. Such a state has minimal 
uncertainty of the spin components transverse to the spin 
quantization axis UM. An explicit construction for the 
coherent state reads |6|: 



|n) = ex.p(iJ z ip) exp(i J x 9) exp(— i J z <j>)\z), 



(8) 



where 8 and (j> are spherical coordinates of n; \z) is the co- 
herent state with the direction of quantization axis along 
z-axis. This definition assures single-valuedness of the 
spin wavefunction. An adiabatic motion of classical spin 
n(t) can be described by the coherent state |n(i)) ac- 
companied with a phase factor e* 7 of purely geometrical 
origin Namely, if the spin moves adiabatically along 
any path I on the unit sphere S 2 of n, the geometrical 
phase 7(7) is equal to a linear integral: 
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(9) 



The local change of the phase is described by Berry's 
connection = (n|id/<9a; M |n). This vector field has two 
components on S 2 : 



(10) 



i a (i-cos( 

-— n i— n) = J — - 

sm dip sm 



The connection A, as well as the geometric phase, is not 
gauge invariant. At a local gauge transformation |n) — > 
exp(zA(n))|n) they are transformed as follows: A — > A + 
dA, 7 — > 7+A/ — Ai , where A is an arbitrary differentiable 
function on S 2 , i and / are its values at the initial and 
final points of path I respectively. However, the phase 
becomes gauge-invariant if the path is closed: I = dc, 
where c is a surface supported by I. In this case: 



7(c) 



dc 



dA = J sin((9)d(9d<j!> = Jfl , (11) 



where f2 is the solid angle subtended by dc at the origin 
of the unit sphere. The integrand in ([ll]) is the field- 
strength B = Jr/r 2 . This field is identical to magnetic 
field produced on S 2 by Dirac's magnetic monopole with 
the charge J located in the center of the sphere. Thus, 
following Berry, we formulated the problem of a localized 
large spin in terms of a scalar charged particle confined 
on the sphere in the field of magnetic monopole. 

In the presence of a crystal field a further simplifica- 
tion becomes possible. As it was shown in Section ||, it 
leads to the localization of the low-energy states near the 
"easy" directions or minima of the field and lowers the 



dimensionality from 2J+ 1 to N, where N is the number 
of the easy positions. The spin trapped near one of the 
easy directions can tunnel to the neighboring minima. 
The tunneling trajectories are solutions of the classical 
equations of motion with imaginary time or velocity. The 
amplitude Wij for the tunneling from the state \i) to a 
neighboring state |j) can be written as my = w expiry. 
Here if is a real, exponentially small factor (see its cal- 
culation in Section ||) and 0y is the Berry's phase along 
the tunneling trajectory connecting the points i and j. 

The set of Berry's phases <^>,j along the tunneling tra- 
jectories {i,j} connecting extrema labeled by i and j 
must satisfy a set of equations. Namely, let us consider a 
plaquette c on the sphere bounded by k tunneling paths 
{h, {i 2 , h}—{ik, h}- Then: 



7(c) = J(fi(c) (mod 4tt)) (12) 



where ik+i — ii and fi(c) is the solid angle subtended 
by the contour dc. The system (|lj) is extended over all 
independent plaquettes. Without loss of generality it is 
possible to consider equations ( |l2| ) only for minimal (ele- 
mentary) plaquettes, i.e., plaquettes of the minimal non- 
zero area whose boundaries do not have self-intersections. 
Equations ([l2]) do not define the phases <fiij unambigu- 
ously. There remains a freedom of a discrete gauge trans- 



formation « 



;+/,— fj containing N real parameters 



fi . One of them can be treated as a common phase fac- 
tor and is inessential. The Schrodinger equation in this 
representation reads: 



n\ii>) = £W 



(13) 



where = YljLi c j\j) ^ s a vector in the iV-dimensional 
space spanned onto the basis \j), j — 1,2...N and H is 
an N x N matrix whose diagonal components are equal 
to a single-well energy level and nondiagonal elements 



are {H} t 



Further we put the diagonal matrix 



elements of H. to be zero. Then equation ( |13| ) can be 
rewritten in the vector form: 



N 

E 



Win Cj 



we* 



(14) 



Eqn. ( |l4| ) is obviously invariant with respect to the 
discrete gauge transformation Wij — > tuye*^* - -^; Cj — > 
c i£ ■ Therefore, any set of phases faj satisfying eqns. 
dl2| ) can be used to find the spectrum and the eigenstates. 

We have seen already that the problem of the quasi- 
classical spin is equivalent to the problem of a charged 
particle confined on the sphere S 2 in the homogeneous 
magnetic field of the monopole. It is a direct spheri- 
cal analog to the problem of a charged particle moving 
on a plane in a homogeneous magnetic field, perpendic- 
ular to the plane. Restricting ourselves with the local- 
ized states, we consider a problem which planar analog 
is the problem of a charged particle living on a 2-D lat- 
tice placed into homogeneous magnetic field. It is known 
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as Harper equation |lq] . The main difference from this 
famous problem studied by Harper, Azbel, Hofstadter, 
Thouless, Wiegmann and many other authors [15 16 1, is 
that, in our case the lattice is embedded into a sphere 
which is a compact manifold, in contrast to the planar 
case. Nevertheless, many features of the Harper equa- 
tions will be encountered here, e.g., sudden variations in 
spectrum at a transition from a rational to an irrational 
flux through an elementary plaquette. 

The initial Hamiltonian 7Y(J) is assumed to possess a 
point group symmetry. It should be noted that 7i(J) is 
invariant with respect to the inversion transformation: 
J — * —J, whereas the reduced effective Hamiltonian is 
not. The reason is that this invariance which stems from 
the time-reversal symmetry cannot be extended onto the 
quantum permutation relations: [Jj, Jk] = j^jkiJi- The 
time reversal requires also anti-linear transformation of 
the state-vectors which cannot be incorporated into 
linear symmetry group. Thus all groups of transforma- 
tions under study consist of rotations only. The point 
groups in 3-dimensions have been studied thoroughly (see 
for example jnj). A special interest will be paid to the 
following point groups: D n , n = 2, 4, 6, O (octahedron), 
and Y (icosahedron) . 

In the next section we show that the action of the 
symmetry transformations onto the effective Hamiltonian 
is not trivial due to the Berry's phases. 



IV. GROUP THEORY ANALYSIS. 

A. Construction of the main representation. 

Let G, a discrete subgroup of 5*0(3), be the point 
group of the crystal field, i.e., the point group leaving 
function /(J) invariant. It always includes the space in- 
version J as a consequence of the time-reversal symmetry. 
Also, we introduce a subgroup G of the full symmetry 
group G = G x d (Ci — {E, I}) which includes rotation 
elements only. Further we employ the notation "symme- 
try group" namely for G. Each group G has several sets 
of equivalent symmetric directions defined by the inter- 
section of equivalent p-fold symmetry axes with the unit 
sphere. Let us denote such a configuration C(G,p) and 
corresponding number of symmetry directions N{G,p) 
(we denoted it earlier as iV). It can be readily seen that 
N(G,p) = \G\/p, where \G\ is the rank of the group G, 
i.e., the number of its elements. The set of N = N(G,p) 
localized states corresponding to the configuration 
G (G, p) is the vector space for a linear unitary represen- 
tation of the group G. This representation depends also 
on J. Let us call it the main representation and denote 
it W(G,p, J). Its dimensionality is obviously N{G,p). 
For J = 0, W(G,p, J) is a matrix representation of some 
subgroup P of the permutation group Sn. Each element 
g of G can be put in one-to-one correspondence to a per- 
mutation "P(g) 6 P. The Hamiltonian of the system is 



invariant under their action. For J ^ or J not equiva- 
lent to 0, the problem becomes quite peculiar since, due 
to Berry's phase factors, the Hamiltonian is no longer 
invariant under the action of the transformations V(g): 



VHV T =W + n. 



(15) 



Hamiltonian TC differs from TL by a gauge transforma- 
tion. Therefore, it is possible to append such a gauge 
transformation IA € U (unitary diagonal matrix) to each 
rotation that the Hamiltonian remains unchanged: 



UVUV T U^ = H. 



(16) 



Thus, a proper representation W(G,p, J) of the symme- 
try group for large spin J or for the Harper's equation on 
the sphere consists of operators W(g) = U(g)V(g). Since 
multiplication of each W(g) by an arbitrary phase factor 
does not violate Eq. (|l6|), the matrices in W(G,p, J) con- 
stitute a projective representation of G in general, that 
is: 

W(ffi)W(sa) = c(g 1 ,g 2 )W(g 1 g 2 ), Si, 92 € G, (17) 

where c{g\, g 2 ) is a function on G x G with values in U{\) 
(2-dimensional cochain). 

Now, a question arises whether the factor set c{g\,g 2 ) 
is equivalent to the trivial one: c'(g±,g 2 ) — 1 for any 
91, 52 € G, as it is for the case J — 0. By definition, 
two factor sets c and c' are equivalent if there exists a 
function b(g) on G with values in U(l) (1-dimcnsional 
cochain) such that: 



c(ffi,52) = 



K9192) 



c' {9l,92)- 



(18) 



We checked for finite groups G C S'0(3) that equa- 
tion ( |l8| ) with c'(gi,g 2 ) = ±1 is really satisfied. In math- 
ematical language it means, that the cochain c is a co- 
cycle but not a coboundary for a half-integer spin and it 
is a coboundary for an integer spin [ pj[ . Therefore, the 
factor set is non-trivial in general. It is equivalent to the 
multiplicative factors {±1} (isomorphic to Z 2 ) which is 
a consequence of the Dirac quantization: 2 J = n, n G N. 
This structure of the factor set might have been antici- 
pated since the parameter space of an arbitrary spin is 
not just SO(3) but its universal covering group SU{2) 
which can be obtained as a non-trivial extension of the 
former one: 1 -> Z 2 -> SU(2) -> 5*0(3) -> 1. In our 
case, spin in CEF, the proper group of symmetries is G 
extended by Z 2 : 1 — > Z 2 — > G — > G — > 1. Instead of 
dealing with the projective representations of G, one can 
work with the linear representations of G. An explicit 
construction of G will be given later in this section. In 
other language we must consider double- valued represen- 
tations of G 11 1 for half-integer J. 

The representation W(G,p, J) was constructed for a 
particular gauge, however, one can easily find the re- 
quired representation if the Hamiltonian undergoes a 
gauge transformation: 

UWA^ = H'. 
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Then a corrected representation leaves the Hamiltonian 
invariant: 

W =uwi£. 



B. Classification of configurations. 

Generally speaking, the iV-dimensional main represen- 
tation is reducible. To perform the reduction of the main 
representation we need to find its characters. They are 
found explicitly in Appendix |A|. Here we issue final re- 
sults. For W(G,p,J), elements with non-zero charac- 
ters are: identity E, the rotation through an angle of 2tt 
about an arbitrary axis Q and rotations C| about the 
p-fold axes. 

X (E) = N; X (Q) = N(-1) 2J ; x (Cg) = 2cos(^); 

X (qQ) = 2(-l) 2J cos(^); (q = l ) ... )P -l) (19) 

Now we proceed to consideration of different point 
groups and their configurations of extrema. 

1. Configurations of the octahedron group O 

Here, we classify possible configurations C(0,p) of the 
octahedron symmetry group O. In general, i.e., with- 
out accidental degeneracy, the minima (maxima) of the 
potential are located either on the equivalent symmetry 
axes of the cube or completely away from them (asym- 
metrically): 

C(0,4) three axes of the fourth order passing through 
the centers of opposite faces, N = 6. 

C(0, 3) four axes of the third order passing through op- 
posite corners, N = 8. 

C(0,2) six axes of the second order through the mid- 
points of opposite edges, N = 12. 

C(0,1) none of the symmetry axes passes through the 
minima N = 24 or N = 48. 

The representations W(0,p, J) of the octahedron group 
acting on the spaces of states corresponding to the above 
described configurations are respectively 6-, 8-, 12-, and 
24(48)-dimensional. 

For a configuration C(0, p) only elements have non- 
zero characters which were calculated earlier. They must 
be divided into classes of conjugate elements. The classes 
with non-zero characters (except E and Q) are: six ro- 
tations C4 and Cf, and three rotations Cf for C(0,4); 
eight rotations C3 and Cf for C(0,3); six rotations C2 
for C(0, 2); and none for C(0, 1). 

The characters are periodic functions of J with the 
period equal to p. It means that the multiplicities of the 



eigenvalues of the Hamiltonian have the same periodicity. 
The characters are invariant under the transformation 
J — > — J (mod p) (reflection). 

The irreducible components contained in representa- 
tions Wjv, N = 6, 8, 12, 24 of the octahedron group are 
given in Table | for values of J inequivalent under the 
translations over p and the reflection. For simplification 
we denoted W(0,p, J) as W N , where N = N(0,p). The 
irreducible components of W^s are not listed since there 
are twice as many of them as those for This rela- 

tionship is correct for representation W\q\ of any group 
G. The characters of the accidental configurations, such 
as the 14-fold configurations on the boundary between 
C(0,4) and C(0,3), are merely sums of the characters 
of the constituting components and can be found from 
the given tables for the basic configurations. 



TABLE I. Irreducible components of the cubic representa- 
tions W N , N = 6, 8, 12, 24. 



J 


C(0,4) 


C(0,3) 


C(0,2) 


C(0,1) 





Ai, E, Fx 


Ax, A 2 
Fi, F 2 


A u E 
Fi, F 2 {2) 


Ai, At, E{2) 
Fi(3), F a (3) 


1 


Fi, F 2 


E, Fi, F 2 


A 2 , E 
Fx (2), F 2 




2 


As, E, F 2 








1/2 


E[,G' 


E[, E' 2 , G' 


E[,E' 2 , G'(2) 


Fi(2), E' 2 (2) 
G'(4) 


3/2 


E' 2 ,G' 


G'(2) 





G 



2. Configurations of the icosahedron group Y 



3. Configurations of Y) 2 



The classification of configurations C(Y,p) for the 
icosahedron group of symmetries is similar to that of the 
octahedron group. Extrema can be located either along 
the directions of the symmetry axes or asymmetrically: 

C(Y, 5) six axes of the fifth order passing through oppo- 
site corners of the icosahedron, N = 12. 

C(Y, 3) ten axes of the third order passing through the 
centers of opposite faces N = 20. 

C(Y, 2) fifteen axes of the second order through the mid- 
points of opposite edges N — 30. 

C(Y, 1) none of the symmetry axes passes through the 
minima N = 60 or 120. 

The main representations of the icosahedron group acting 
on the spaces of the configurations are respectively 12-, 
20-, 30-, and 60(120)-dimensional. The classes with non- 
zero characters, besides E and Q, are: twelve rotations 
C5' 4 and twelve rotations C5' 3 for C(Y,5), twenty rota- 
tions C 3 ' for C(Y,3), fifteen rotations C2 for C(Y,2), 
and none for C(Y, 1). The multiplicities of the eigenval- 
ues of the Hamiltonian, for a configuration C(Y,p), have 
the period p. The irreducible components contained in 
representations Wn, N — 12, 20, 30, 60 of the icosahe- 
dron group are given in Table [fi| 



Next, we consider the configurations of three groups of 
symmetries Dat (N — 2, 4, 6). Despite their simplicity, 
Berry's phase introduces here some interesting effects as 
well. 

The configurations of D2 are quite simple: 

C(D2, 2) one axis of the second order, N = 2. 

C(T>2: 1) none of the symmetry axes passes through the 
extrema, N — 4 or 8. 

The characters of the D2 representations and the irre- 
ducible components contained in representations Wjv, 
TV = 2, 4 are given in Table III. 



TABLE III. 


Irreducible components 


of the D2 representa- 


tions. 






J 


C(D 2 ,2) 


C(D 2 ,1) 





A, B 3 


A, Bi, B2, B3 


1 


B%, B 2 




1/2 


E' 


E'(2) 



TABLE II. Irreducible components of the icosahedron rep- 
resentations W N , N = 12, 20, 30, 60. 



J 


C(Y,5) 


C(Y,3) 


C(Y,2) 


C(Y,1) 




1 

2 


A, Fi 
F a , H 

Fx, G, H 

Fi, G, H 


A, Fx, F 2 
G(2), H 
Fx, F2 
G, H{2) 


A, Fx, F 2 
G{2), H(3) 
Fi(2), Fa (2) 
G(2), H{2) 


A, Fx (3), F 2 (3) 
G(4), H(5) 


1/2 

3/2 
5/2 


E[, G', I' 
E' 2 , G', I' 

I'm 


E'x, E' 2 
G', I'(2) 
G'(2), J'(2) 


E[, E' 2 
G'(2), 


E'x{2), E' 2 {2) 
G'(4), J' (6) 
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4- Configurations of the tetragonal group D4 



5. Configurations of the hexagonal group De 



The configurations of D4 are listed below: 

C(T>i, 4) one axis of the fourth order, N = 2. 

C(D4, 2) two axes of the second order, N = 4. 

C(D4, 1) none of the symmetry axes passes through the 
extrema, N = 8 or 16. 

The irreducible components contained in representations 
W N , N = 2, 4, 8 of D 4 are given in Table |iy|. 

An interesting conclusion can be drawn from the data 
in Table [V. The 2- fold classical degeneracy of the con- 
figurations of C(D4, 4) is not lifted for all but even values 
of J. Thus, the tunneling is allowed only for even spins. 
This result cannot be accounted for by Kramers degener- 
acy, as it was possible in Q for D 2 configuration, and is 
totally due to the symmetry combined with the Berry's 
phase. Also, it shows importance of the details of the 
background, i.e., D2 (considered in the previous section 
and in iQ), D 4 , and D 6 (considered in the next section) 
groups of symmetries have the easy axis (2-fold) con- 
figuration, however the tunneling is allowed in the Djy 
environment only for J = (mod N/2) and is defined 
by the anisotropy in the plane normal to the easy axis. 



Due to similarity of this group with D4, we just present 
the data on the T)$ representations. 

C(Dq, 6) one axis of the sixth order, N = 2. 

C(Dq, 2) three axes of the second order, N = 6. 

C(Dq, 1) none of the symmetry axes passes through the 
minima, N = 12 or 24. 



TABLE V. Irreducible components of the D6 representa- 
tions. 



J 


C(D 6 ,6) 


C(D 6 ,2) 


C(D 6 ,1) 





Ai, A 2 


A\, Bi, Ei, E2 


Ai, A 2 , Bi, B 2 , Ei(2), E 2 {2) 


1 


Ej 


A2, B2, Ei, E2 




2 


E 2 






3 


B\, B2 






1/2 


E[ 


E[, E'2, E'z 


E[{2), E l 2 {2) , E' z (2) 


3/2 


E' 3 






5/2 


E' 2 







TABLE IV. Irreducible components of the D4 representa- 
tions. 



J 


C(D 4 ,4) 


C(D 4 ,2) 


C(D 4 ,1) 





Ai, A 2 


Ai, Bi, E 


Ai, A 2 , Bi, B 2 , E(2) 


1 


E 


A 2 , B 2 , E 




2 


Bi, B 2 






1/2 


E[ 


E'i, E' 2 


E'i{2), E' 2 {2) 


3/2 


E' 2 
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V. SPECTRUM 

The group-theoretical analysis of the last section gives 
the number of splitted sublevels in the initial ./V-fold mul- 
tiplet and their degeneracies. In this section we find the 
order of the sublevels and distances between them. It 
requires explicit diagonalization of the reduced Hamilto- 
nian. As we show below, the spectrum is much more 
subtle matter than the number and degeneracy of the 
sublevels. It may depend on details of the Hamiltonian. 

We assume that all tunneling paths between nearest 
minima are equivalent, that is, all non-zero tunneling am- 
plitudes have equal absolute values \w\. Consequently, 
w enters the Hamiltonian as a common multiplier and 
all eigenvalues are multiples of w in zero magnetic field. 
The solid angle covered by the minimal non-trivial closed 
path will be assumed known. It is, actually, a constant 
for all configurations but C(G, 2), G = O, Y, where it 
is a function of some dimensionless combinations of the 
CEF parameters, e.g., ratio u in Section 0. 

In some cases, not only in simple ones, such as 
C(T>n,N), there are two tunneling trajectories connect- 
ing nearest minima. E.g., in a vicinity of the boundary 
between 6- and 8-fold configurations of the cubic group 
(see Fig. ||), the tunneling trajectory deviates from the 
geodesies connecting the minima and, due to the sym- 
metry, there are two trajectories located, symmetrically 
with respect to the geodesies. However, the two trajecto- 
ries can be considered as one effective path with the tun- 
neling amplitude of 2w cos( JO/2) (see Eq. (pcf)), where 
w is the tunneling amplitude of a single path and O is 
the solid angle subtended by the two trajectories. 

Before proceeding to a detailed analysis of the spectra, 
we obtain some relations between eigenvalues of the same 
configuration, but for different J. These relations are of 
purely geometric origin Let us assume that the pa- 
rameter space S 2 can be covered completely and without 
overlap by s congruent plaquettes whose boundaries are 
the tunneling trajectories. E.g., these are 2 hemispheres 
for C(Djy,2), N = 2, 4, 6, configurations; N orange-like 
segments for C(Djv, N), Fig. [|; 8 curved right-angled tri- 
angles for C (0 , 4) , Fig. H. Each plaquette subtends a solid 
angle of Att/s, and Berry's phase for each loop is AirJ/s. 
Then, from fll2]), it follows that the spectrum is a peri- 
odic function of J with the period s/2 FT. The spectra of 
systems differing by transformation 7(c) — > —7(c) must 
be identical due to the time-reversal symmetry. Hence, 
all J's are divided into s/2 + 1 equivalence classes defined 



x This statement is conventional: it is periodic if w does not 
depend on J. However, the ratios of the interlevel distances 
are periodic functions of J. 



by a set of numbers 0, 1/2, 1, . . . , s/4. A fixed J belongs 
to the class of equivalence labeled by: 

min|J + ns/2|. (20) 

Hereinafter, we will work only with the minimal non- 
equivalent J's. 

In a more general setting, i.e., in the presence of n 
different elementary plaquettes, periodicity of the spec- 
tra depends on the rationality of the flux quanta passing 
through each plaquette: if a flux per each plaquette is 

= Jfli — 2-nJPi/Qi, i — 1, . . . , n, where Pi and Qi are 
mutually prime integers, then the period of the spectra 
is the least common multiple of Qi, i = 1, . . . , n. Other- 
wise the spectra are not periodic and each J represents a 
class. Thus, if n = 1 the spectra is always periodic and if 
n > 1 it is not in general (unless an additional symmetry 
is present). 

An extra symmetry of the spectra can be extracted by 
considering an operation of the change of sign: w — ► — w. 
This transformation inverts energy levels inside of each 
class. On the other hand, the spectra depend only on 
gauge invariants w k cos(Jttk), where ttk = Amir/s (m is 
an integer) is the solid angle subtended by a closed con- 
tour containing k tunneling paths and is a multiple of the 
solid angle subtended by the elementary plaquette Att/s. 
If all closed contours contain even number of the paths (k 
is even), e.g., C(0, 3), C(D 4 , 2), the levels are symmetric 
inside of each class, that is, they come in pairs of oppo- 
site sign ±E. For example, for eigenvalues of C(0,3) 
the following relations are satisfied: E(J — 0;Ai) = 
-E(J = 0;A 2 ), E(J = 0,1; Fx) = -E(J = 0,1; F 2 ), 
E(J = 1;E) = 0, E(J = 1/2; E[) = ~E(J = 1/2; 7^), 
E(J = 1/2; G") = 0. If some of the closed contours con- 
sist of odd number of the paths, e.g., C(0,4), C(Y,5), 
then the simultaneous change of sign w — > — w and 
shift J — ► J + s/4 leaves the invariant combinations un- 
changed. Therefore, each level E in the class of J has 
its counterpart — E in the class of J + s/4. For ex- 
ample, in C(0,4): E(J = 0;Ai) = -E(J = 2;A 2 ), 
E(J = 0;^) = -E(J = 2;F 2 ), E(J = 0;E) = 
~E{J = 2;E), E(J = 1/2; E[) = -E(J = 3/2; E' 2 ), 
E(J= 1/2; G ) = -E(J = 3/2; G'). If J and J + s/4 be- 
long to the same equivalence class, their spectrum is sym- 
metric, e.g., in C(0,4): E(J = l;Fi) = -E(J = 1;F 2 ). 

Further in this section we calculate the spectra for dif- 
ferent groups of symmetry and configurations. 

A. Spectra of the D„ (n = 2, 4, 6). 

The configurations of D4 are shown on Fig. |[ In the 
case of the Dg configurations, there are six minima on the 
equatorial circle (C(Dq, 2)) and six tunneling paths con- 
necting the antipodal points (C(D 6 ,6)). For C(D 2 ,2), it 
is just two minima connected by two tunneling trajecto- 
ries. The total tunneling amplitude for C(Dn, N), from 
one pole to the other is: 







iV-1 



■w 



^2 exp(iA-KkJ/N), 



(21) 



k=0 



where we prescribed a phase factor of unity to one of the 
tunneling paths. The Hamiltonian is a 2 x 2 matrix with 
the following eigenvalues: 

± 2wcos(ttJ) for N = 2, 
E = { ± 4ujcos(ttJ) cos(ttJ/2) for N = 4, 

± 2u;cos(ttJ)(1 + 2 cos(2tt J/3)) for N = 6. 

(22) 

In full agreement with the predictions of Section the 
paths interfere destructively for all spin values but J = 
(mod N/2). 

The case of minimal symmetry C(D2, 2) has been con- 
sidered by D. Loss et al (T^j earlier. They argued that 
in the case of half-integer J the tunneling amplitudes 
along the two paths cancel each other. One can see from 
Eqns. ( p2| ) that, when the number of equivalent tunnel- 
ing paths increases due to the symmetry, such a cance- 
lation takes place for integer J as well (with exception 
of J = (mod N/2)), where the classical degeneracy of 
the ground-state level is 2-fold for all C(Djv, N). 

In the presence of magnetic field the eigenvalues are 
E(h) = ± V / E 2 (0) + (hJ) 2 , where h = gfx B H and H is 
the component of magnetic field along the easy direction. 

For the C(Djv,2) N = 4, 6 configurations, the Hamil- 
tonian is that of the one-dimensional iV-site tight-binding 
model (|] , with eigenvalues 0: 

E k = 2wcos(2n(k + J)/N), k = 0, 1, . . . , N - 1. (23) 

The magnetic field enters the Hamiltonian as a site- 
diagonal matrix: 

H h = -hJcos{<j) h - 2-kI/N), I = 0, 1, . . . ,N - 1, (24) 

where h = gfi^H , H is the in-plane component of mag- 
netic field, and 4>h is the angle of this component with 
respect to the easy direction of the CEF labeled by / = 0. 
The eigenvalues oiTl + Hh can be found analytically. For 
C(D4, 2), one finds: 



classes without a magnetic moment in the ground-state 
(integer spins) and has a Curie-like behavior for the ones 
with a magnetic moment in the ground-state (half-integer 
spins); see Appendix |c| for details. 

In the case of C(T>n, 2) configurations, there is a spec- 
tral difference between integer and half-integer spins only 
which can be ascribed to Kramer's degeneracy. In the 
next section, we consider non-Abelian cases, where more 
complex division on equivalence classes occurs. 



TABLE VI. Spectrum of C(D 4 ,2) in magnetic field, the 
limit of small magnetic field, and low temperature magnetic 
susceptibility (J3 = l/(k B T). 



J 


Eigenvalues 


Susceptibility 




1/2 


±(2w + \h z J' 2 /w), ±\h A J' A wa.(2<t> h )/w 
±(V2w ± \hJ + ^h 2 J 2 /w) 


^{gJ^f/w 



TABLE VII. Spectrum of C(D 6 ,2) in magnetic field, the 
limit of small magnetic field, and low temperature magnetic 
susceptibility. 



J 


Eigenvalues 


Susceptibility 





±{2w + J' d /w) 
±(w-'^h 2 J 2 /w) 
±{w + lh 2 J 2 /w) 


(gj/isf/w 


1/2 


(2), ±(v / 3w ± \h.J + ^§h 2 J 2 /w) 


±(gJgB) 2 /(k B T) 



E z = 2w 2 



~±]J^cos 2 (TTj) + 2h 2 w 2 + tjcoa 2 (2<f> h ) 



where we used h as a shorthand for hJ. The spectra 
of C(D/i, 2) (previously calculated in ||) and C(Dq, 2) in 
magnetic field are giv en in the limit of small magnetic 
fields in TablesjvJ and VH respectively. The last column 
of Tables VI, VII is the low temperature magnetic sus- 
ceptibility which is a readily observable physical quan- 
tity. The susceptibility saturates to a constant for the 



^The label k in Eqn. 
value of the level. 



j does not correspond to the algebraic 
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1/2 
3/2 



-V2w(4), 2 v / 2w(2) 
-2v / 2w(2), \/2w(4) 



2/3/9 
0/9 




FIG. 4. The configurations of D4. Filled circles and dashed 
lines belong to C(D4,4), and filled squares and solid lines 
belong to C(D4,2) respectively. 



B. Spectra of the O configurations 

The cubic symmetries are quite common in nature. We 
will perform a detailed study of the configurations of the 
octahedron group. The Hamiltonian for configuration 
C(0,4) has the following matrix elements: 

hu = 0, i = l, 2,... 6, 

hij=0, \i-j\=l,i + j = 3, 7,11, (25) 
I I = |w|, for other 1 < i, j < 6 

where we adopted the enumeration shown on Fig. |l|. The 
tunneling trajectories divide the sphere into eight plaque- 
ttes. Relation (|l2|), written for each plaquette, gives eight 
equations for the phases 0y ;, where 0(c) = tt/2 (an exam- 
ple of a set of the phases for this configuration as well as a 
calculation of the spectra is given in Appendix |b|) . Only 
seven equations are independent. Given definite phases, 
the diagonalization is straightforward. The eigenvalues 
can be expressed in the following closed form p3]: 

E k (J) = {-l) k 2w X (n(J + 2k)), k = 0, . . . , 5, (26) 



y(x) = cos — cos cos — h sin — sin 

AW 3 2\ 3 3 2 



1/2 



The ordered spectra of C(0,4) are given in Table VIII 
(w > 0) for the minimal set of J's; the spectra for other 
J's can be obtained by the use of the equivalence rela- 



tion ( 20 ) . Note that the spectra should be inverted if w 



is negative. 

TABLE VIII. The spectra and the low temperature mag- 
netic susceptibilities of C(0,4) (/3 = l/(ksT), a common fac- 
tor of (<7J/xb) 2 is omitted ). 

J Eigenvalues (degeneracies) Susceptibility 

() -2w (2), 0(3), Aw (I) l/(3u>) 

1 -2w (3), 2w (3) /3/6 

2 -4w(l), 0(3), 2w(2) l/(6w) 
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The physical difference among the classes is manifested 
when magnetic held is applied. The magnetic part of the 
Hamiltonian in this case is: 

Hh = J diag(-h z , h z , -h x , h x , -h yi h y ). 

The full Hamiltonian can be easily diagonalized for some 
symmetric direction of the held, e.g., along easy direc- 
tion (1,0,0). The direction of the held does not influ- 
ence the low temperature susceptibility since the latter 
is isotropic in a cubic CEF. However, individual levels 
of the ground-state multiplet may have anisotropic mag- 
netic susceptibility as well as anisotropic magnetization. 
The low temperature magnetic susceptib ilities of C(0, 4) 
are collected in the last column of Table VIII . 



J 


Eigenvalues (degeneracies) 


Susceptibility 





-3w(l), -w(3), w(3), 3w(l) 


l/(3u>) 


1 


-2w (3), 0(2), 2w (3) 


0/6 


1/2 


-V6w (2), 0(4), v%™(2) 


0/9 


3/2 


-V3w(4), y/3w(A) 


20/9 




FIG. 5. The paths of the spin on the unit sphere between 
the easy positions of the field. The case of C(0,3) configura- 
tion. 

In the case of C(0,3) configuration, the minima are 
located at the ver tices of a cube inscribed into the unit 
sphere: sind = -\/2/3, sin(2</>) — 0, where 9 and (f> are 
the spherical coordinates of the minima (see Fig. |). The 
tunneling trajectories divide the surface of the sphere into 
six congruent plaquettes; each subtends a solid angle of 
27r/3. Five independent equations ( |l2| ) hx the tunnel- 
ing phase shifts and the Hamiltonian up to an arbitrary 
gauge transformation. The eight eigenvalues are pM: 



E± = ±2<(tt( J + 3k)), k = 0, 1, 2, 3, 



3 + 2 cos x cos h 4 cos — cos — g(x) 

3 2 3 w 



g(x) = (4sin 



2 ^ -2 

- sir 



1)< 



(27) 



The ordered eigenvalues are presented in Table [X for 
the non-equivalent J's (w > 0). Analysis of the mag- 
netic response is quite straightforward as well (see Ap- 
pendix |Cj); the magnetic susceptibilities of the classes are 
given in the last column of Table IX . 



TABLE IX. The spectra and the low temperature magnetic 
susceptibilities of C(0,3) (0 = l/(k^,T), a common factor of 
(gJg,B) 2 is omitted). 
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Consideration of C(0,2) will be postponed till Sec- 



tion VI 



1. Multiple tunneling path regime. 

In C(0,4) configuration, a tunneling trajectory con- 
necting two minima, e.g., minima 3 and 5 on Fig. 0, is 
not necessarily a geodesies on the sphere. For example, if 
the mid-point of the geodesies connecting minima 3 and 
5 is a maximum of the CEF potential then the tunneling 
trajectory connecting the minima will split in two paths: 
one deviating towards the "north" pole (minimum 1) and 
the other towards the "south" pole (minimum 2) as it is 
shown on Fig. |^. One path is a mirror copy of the other 
with respect to the "equatorial" plane. Thus, the ab- 
solute values of the tunneling amplitudes corresponding 
to the two trajectories (\w\) are identical. To find the 
compound tunneling amplitude we assume that one of 
the trajectories, e.g., the one connecting minima 3 and 5, 
and located in the "south" hemisphere, has the phase (p%: 
wi = \w\ exp(i(pi). Then, due to the Berry connection, 
the other amplitude must be w 2 = \w\exp{i{ip\ — J£t)), 
where ft is the solid angle subtended by the two trajec- 
tories. The effective amplitude is: 

w e = w x +w 2 = 2\w\e i ^ 1 - jn W cos(jn/2). (28) 



Interesting conclusions can be derived from formula (28). 
Firstly, the splitting of the trajectories does not change 
the connectivity matrix of the configuration it just mod- 
ifies the multiplier of Hamiltonian (|25|) and all results 
obtained for configuration C(0,4) hold true. Secondly, 
the spectrum may be an oscillating function of J or, if 
one would be able to vary parameters in such a way that 
n changes from its maximum value to zero, several os- 
cillations of the spectrum could be observed as well. To 
estimate the number of oscillations we use the fact that: 
different tunneling trajectories emanated from a site and 
ending at some other site(s) do not intersect at interme- 
diate points (they can only intersect at the end points). 
Then we can state that the maximal possible deviation of 
the trajectories from the spherical geodesies connecting 
the positions of C(0, 4) configuration is reached when the 
trajectories pass along the spherical geodesies connecting 
the geometrically closest positions of C(0, 3) and C(0, 4) 
configurations. Fig. [?] depicts this situation: the two tun- 
neling trajectories connecting the 6-fold global minima 3 
and 5 (filled circles) are passing very closely to the 8- 
fold local minima (filled triangles), thus, "avoiding" the 
12-fold global maxima (filled hexagons). The solid angle 
enclosed by the two trajectories (shaded area on Fig. [?]) 
varies in the range < 57 < 7r/3. Upon such a variation 
of fl the spectrum will make J/12 full oscillations. 




FIG. 6. Tunneling trajectories of the spin (double paths). 
The 6- fold configurations of O. The hexagons show the loca- 
tions of the maxima of the CEF potential. 




• - 6-fold (global minima) 
▼ - 8-fold (local minima) 
#- 12-fold (global maxima) 

FIG. 7. Tunneling trajectories of the spin (double paths). 
The 6- fold configurations of O. The tunneling trajectories 
pass closely to the local minima (locations of the 8-fold con- 
figuration) . 

Next we analyze the multiple tunneling trajectories 
of C(0,3). Fig. H depicts the splitting of the trajec- 
tory connecting minima 1 and 5 (solid curves A and B). 
The situation is similar to that of C(0,4) configuration 
(the oscillations take place and their maximal number is 
J/12) except one subtle point: when a trajectory devi- 
ates strongly from the geodesies it approaches the trajec- 
tory connecting a next-nearest-neighbor (dashed lines on 
Fig. ||), e.g., lines A' and B' which connect 1 with 4 and 
8 respectively. This is a very drastic change in the tun- 
neling regime which leads to a change of the connectivity 
matrix. 
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FIG. 8. Tunneling trajectories of the spin (double paths). 
The 8- fold configurations of O. The hexagons show the loca- 
tions of the maxima of the CEF potential. 

To calculate the spectrum we assume that the ab- 
solute values of the single tunneling amplitudes to the 
nearest- and next-nearest-neighbor sites are the same w. 
However, the effective amplitude for the nearest-neighbor 
tunneling is 2w cos(Jfl/2) due to the double trajectories. 
The elementary plaquette, in this case, is a triangle cov- 
ering the solid angle of 7r/3, e.g, triangle 1-5-8-1 on Fig. ||. 
In this case, plaquettes cover the sphere twice. Then, the 
periodicity of the spectra is given by s = 4ir/(ir/3), which 
is a half of the total number of the elementary plaquettes. 
Application of the symmetry arguments given at the be- 
ginning of this section leads to the following properties of 
the spectrum: the periodicity of the spectral behavior is 
J = s/2 = 6, |J + 6n|, n £ Z is equivalent to J, the spec- 
trum of J + 3 is the inverted spectrum of J. The results 
of the diagonalization are summarized in Tabic [x|. 

TABLE X. The spectra of C(0, 3), the region of the mul- 
tiple tunneling path regime [x = 2cos(JT2/2)]; all eigenvalues 
are multiples of w. 

J Eigenvalues (degeneracies) 

1 
2 
3 

1/2 
3/2 
5/2 



-3(1 - x) (1), -(1 + x) (3), (1 - x) (3), 3(1 + x) (1) 

-2(1 - i/2) (3), -3x (2), 2(1 + s/2) (3) 

-2(l + a;/2) (3), 3x(2), 2(1 -x/2) (3) 

-3(1 + x) (1), -(1 - x) (3), (1 + x) (3), 3(1 - x) (1) 

-( V6 - aV3) (2), - xV3(4), ( Ve + xVs) (2) 

-V3(l+a*)(4), ^3(1+^(4) 

-{Vs + xVi) (2), aV3(4), (V6-xVS) (2) 



C. Spectra of the Y configurations 

The analysis of the configurations of Y group is te- 
dious, though similar to that for O group. We present 
only the results of the analysis here. Table |Xl| contains 
the spectra and the low temperature susceptibilities of 
C(Y, 5) configura tion (the energies are multiples of w). 
Tables XII, XIII contain the spectra and the low tem- 



perature susceptibilities of C(Y, 3) configuration respec- 
tively. 



TABLE XL The spectra and the low temperature magnetic 
susceptibilities of C(Y,5); all eigenvalues are multiples of w 
and all susceptibilities are multiples of (gJfiB) 2 [(3 = l/(ksT), 
Ci = cos(7r/10), and C3 = cos(37r/10)]. 



J Eigenvalues(degeneracies) 



-V5(3), -1(5), V5(3), 5(1) 
-V5 (4), (V5 - 3)/2 (5), (5 + v / 5)/2 (3) 
-VE(4), (V5-5)/2(3), (V5 + 3)/2(5) 
-(VE+ 3)/2 (5), (5 - v / 5)/2 (3), ^5 (4) 
-(5 + >/5)/2(3), (3-V5)/2(5), ^5(4) 
-5(1), -V5(3), 1(5), >/5(3) 
-2ci (6), (3 - v/5)ci (4), 2V5d (2) 
-2VEc 3 (2), -2c 3 (6), (3 + VE)c 3 (4) 
-V5(6), V5(6) 

-(3 + V5)c 3 (4), 2c 3 (6), 2V5c 3 (2), 
-2VEa (2), (V5 - 3)ci (4), 2 Cl (6) 



Susceptibility 

(1 + y/5)/(Gw) 

0/9 

(3/9 

2(3/9 

(5 + VE)/(30w) 
(3/5 

(5 + \/5)ci/(15 
(3(5 + VE)/30 
(3/5 
(3/9 



TABLE XII. The spectra of C(Y,3); all eigenvalues are 
multiples of w. 



J 



Eigenvalues (degeneracies 




1 

2 
3 

1/2 
3/2 
5/2 



V5(3), -2(4), 0(4), 1(5), V5 (3), 3(1) 



-(1 + VlS)/2 (5), -1 (4), (3 - v/5)/2 (3), 
(-l + Vl3)/2(5), (3 + V5)/2(3) 
-(3 + V5)/2 (3), (1 - 7T3)/2 (5), (3 - y/§)/2 (3), 
1(4), (l + Vi3)/2(5) 

-3(1), -^5(3), -1(5), 0(4), 2(4), >/g(3) 

-(%/3 + V7)/2 (6), >/3(l - V5)/2 (2), 

(- V3 + v / 7)/2 (6), >/3 (4), >/3(l + y/S)/2 (2) 

->/6(4), -1(6), 1(6), V6 (4) 

-V3(l + v / 5)/2 (2), -V3 (4), (V3 - v / 7)/2 (6) 

v^(-l + 75)/2(2), (V3 + V7)/2(6) 
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TABLE XIII. Low Temperature magnetic susceptibili- 
ties of C(Y,3) [a common factor of {gj^s) 2 is omitted, 
= l/(kaT)]. 



J 




1 

2 
3 

1/2 

3/2 
5/2 



Susceptibility 
(7/6 + 11^5/18) /it; 



P(12y/E + 37 + ^13(3^5 + 4))/468 
_ 0/6 
(V5 + 3)/(6tu) 
0(10^ + 83 + ^27(5^5 - 2))/630 
0(4 v / 5 + 9)/9O 
0/9 



The multiple tunneling path regime is present in con- 
figurations C(Y, 5) and C(Y, 3) as well. Its analysis is 
similar to that of configurations C(0, 4) and C(0, 3). We 
present here its summary only: The regions of existence 
of configurations C(Y, 5) and C(Y, 3) in the parameter 
space of the CEF are divided into two parts for each con- 
figuration. One part corresponds to the single tunneling 
path regime. The above theory is valid in this region. 
The other part is of the multiple tunneling path regime. 
The spectra are oscillating functions of J in this region 
since w ~ cos(Jfi/2), < ft < 27r/15. Upon full mono- 
tonic variation of Q, the spectra makes ~ J/30 oscilla- 
tions for the both configurations. The spectra of C(Y, 5) 
given in Table XI holds valid for the both regimes. For 
C(Y, 3) configuration, in a range of parameters the prox- 
imity of the minima positions may be altered: each min- 
imum position (vertex of the dodecahedron where some 
three faces intersect) should be geometrically connected 
not just to the three nearest-neighbors but also to the six 
next-nearest-neighbors. 



VI. RANDOM ENERGY LEVELS. 



For all configurations considered in previous section, 
the spectra were simple periodic functions of J, which 
was due to the fact that a rational number of flux quanta 
(<f> = Jfi(c) = 2irJP/Q) passes through each plaque- 
tte. This is not the case for more complex configurations 
such as C(G, 2), G — O, Y. In Fig. || we present the 
spatial distribution of minima of C(0,2) configuration. 
The segments connecting the minima are not real tun- 
neling trajectories but rather guidelines. The tunneling 
paths may deviate strongly from the geodesies connect- 
ing corresponding minima both to the locations of the 
6-fold (centers of the cube faces) and 8-fold (vertices of 
the cube) configurations' positions. The exact form of 
the paths depends on the CEF constants, e.g., for the 
simplest Hamiltonian (^|), where configuration C(0, 2) is 
realized, it is a function of ratio b/a. Instead of study- 
ing non-universal tunneling trajectories, we introduce a 
parameter a (0 < a < 2ir/3): the solid angle subtended 
by a square-like contour. The solid angle subtended by 
a triangle-like circuit is 7r/2 — 3a/4. A knowledge of this 
parameter together with w is sufficient to define the spec- 
tra of the 12-fold configuration. Since a may be an ir- 
rational multiple of ir, the spectra as a function of J is 
not expected to be a finite set of values, but a fractal 
set. T he sp ectra of the 12-fold configuration are given in 
Table XTV. The spectra undergo J/12 oscillations upon 
a monotonic variation of < a < 2ir/3 for a given value 
of spin J. 

Configuration C(Y, 2) is even more complex than 
C(0,2). Its minima directions correspond to the mid- 
points of the icosahedron edges (see Fig. [H]). The pa- 
rameter a (0 < a < 7r/3) here corresponds to the solid 



15 



angle subtended by a pentagon-like contour. The spectra 
undergo J/30 oscillations upon a monotonic variation of 
< a < 7r/3. The spectra of the 30 -fold configuration 



for odd values of J are given in Table XV. 

The spectra described in this section have features of 
randomness. Indeed, the function {aJ} (fractional part 
of aJ) with an irrational a is known as a generator of 
random numbers. Thus, the ratios of the transition fre- 
quencies for configurations C(G, 2), G — O, Y vary in 
an uncontrollable way when large J changes by 1. This 
behavior differs dramatically from that for other cubic 
and icosahedral configurations which display permanent 
ratios of the frequencies for a fixed configuration. Thus, 
the configurations C(G, 2), G = O, Y realize the chaotic 
spectra of deterministic systems. This situation is well- 
known, e.g., for the Hydrogen atom in a uniform mag- 
netic field p9| . The peculiarity of our problem is that it 
displays chaos in a finite set of numbers (12 or 30) and 
that the chaotic behavior can be found analytically. An- 
other special feature of our system is that stochasticity 
in it is combined with deterministic multiplicity distri- 
bution. For example, in the case of the C(0,2) configu- 
ration the 12 levels are divided into submultiplets given 
in Table [|, independently on a. However, their mutual 
arrangement is unpredictable. 

For a two parametric Hamiltonian, e.g., Hamilto- 
nian (|J) for the octahedron group, the configurations 
C(G, 2), G = O, Y correspond to the single tunneling 
regime. The multiple tunneling regime may occur if the 
invariants of higher orders are included. 



TABLE XIV. Spectra of C(0, 2). All eigenvalues are mul- 
tiples of w; x = J (a + 27r)/4. 



Integer J 


Half-integer J 


Energy (Degeneracy ) 


Energy (Degeneracy) 


4 cos a; (1) 


2(cosa;± v^sins) (2,2) 


-2cosz (2) 


— cos X ± \/2 + cos 2 x (4,4) 


2 cos x (3) 




— cos x ± \/8 — 7 cos 2 x (3,3) 






FIG. 9. Minima distribution of the 12-fold configurations 
of O. 




FIG. 10. Minima distribution of the 30-fold configurations 
of Y. 



TABLE XV. Spectra of C(Y,2) for odd values of J. All 



eigenvalues are multiples of w, x = cos (J (a 


+ 3tt)/5). 


Energy (Degeneracy) 


±1 + 2x (4,4) 




-a; ± V4 - 3a; 2 (5,5) 




(1 + VS){-x ± V 4 + ( 5 - 4 v / 5)a; 2 


)/2 (3,3) 


(1 - VE)(-x ± V4 + (5 + 4 v / 5)a; 2 


)/2 (3,3) 



1G 



In the presence of infinitely small magnetic field the 
ground-state of configuration C(0,2) acquires either a 
finite magnetic moment or a finite susceptibility. We an- 
alyzed this problem for the field directed along one of the 
fourth order axes and w > 0. Then the finite magnetic 
moment 2gJ^\ sinx|/(2(8 — 7 cos 2 a;)) 1 / 2 is acquired at 
cos a; > —1/2 (x = J (a + 2n)/4), otherwise the finite 
magnetic susceptibility x — — (.9 ^Mb) 2 /(3w cos x) occurs 
for integer J. For half-integer J, the magnetic moment 
<7<//ib/3 is acquired at cosx < cos(37r/8), otherwise this 
value of the moment is multiplied by a factor: 

i 

/ c 2 + 5 + 3CV2T? + 2y/2\ sin(x) | (3c + v / 2+~c 2 ) A 2 



where c = cos(x). 
values. 



2(2 



Note the random character of these 



VII. NUMERICAL ANALYSIS. THE CASE OF 
THE CUBIC CEF. 

The main obstacle to a reliable numerical analysis 
of the problem is the fact that nobody knows how the 
Hamiltonian looks like. The case of the rare-earth ions 
with large total angular momenta interacting with the 
CEF represents an exception. Only the orbital part L of 
the total angular momentum of a single magnetic elec- 
tron interacts with the crystalline field. All terms, in 
the expansion of the crystalline potential with the degree 
larger than 21, where I is the orbital quantum number 
of the single magnetic electron, vanish |2^j , thus, sim- 
plifying the analysis. For the 4/-group electrons with 
I = 3, this gives the highest non- vanishing terms of the 
sixth order. Considering a CEF of a particular symme- 
try group brings further simplification, e.g., in the cubic 
CEF, there are only two independent invariants of the 
sixth order and one of the fours order. The two of the 
sixth order are combined in one invariant (sec Eq. (](])) 
for a real interaction which is the Coulomb interaction 
between the charge carriers. 

It has been shown in Section [il] that Hamiltonian (|^) 
has configurations C(0, 4), C(0, 3), andC(0,2) as sets of 
its classical extrema. In this meaning it is rather general. 
Therefore, we apply numerical analysis to this Hamilto- 
nian in a wide range of J's. It means that we diagonalize 
numerically (2J+1) x (2 J+ 1) matrix for one-parametric 
set of Hamiltonians ( |29|) . The choice of this Hamiltonian 
is partly justified by the above consideration. Our pur- 
pose is to find numerically what J can be considered as 
large, i.e, starting from what J our theory gives satis- 
factory description. The second important problem is 
the crossover behavior of the spectrum near configura- 
tion boundaries described in Section |J. 

First numerical studies of the crystal field effects on 
angular momenta were performed in the early 60's by 



Lea, Leask, and Wolf (20|. These authors studied a cu- 
bic crystal field Hamiltonian similar to ([|). Their main 
interest was how the angular momentum degeneracy of 
/-electrons is lifted. For this purpose it was enough to 
consider values of J spanned from 3 to 8. Refraining 
ourself from this limitations, we study numerically the 
Hamiltonian consisting of terms of the fourth and sixth 
order for an arbitrary value of J. However, we use a 



different parameterization than that used in |20| for the 
same Hamiltonian: 



H? = - 



cos(^)C>2 5 sin(0)Og 



(J(J + 1)) 2 14(J(J + 1))3 



(29) 



where O® and 0§ are Stevens' operator equiva- 
lents pjjpg l, and 4> is a parameter taking values in the 
interval [— tt,tt]. Our parameterization corresponds to 
a unit circle on the phase diagram of Hamiltonian (^) 
(see Fi g. p[ ): a = cos((/>), 6 = sin(0), whereas that cho- 
sen in [pOfl corresponds to the square: a — x, 146/5 = 
±(1 — \x\); — 1 < x < 1. The coefficient of 5/14 reflects 
the difference between our invariant of the sixth order 
in (^]) and the commonly used Stevens' operator equiva- 
lent Og. 

For relatively small values of spins 2 J + 1 ~ N , where 
N is the number of extrema of the CEF, the quasi- 
classical description fails and the spectrum of Hamil- 
tonian ( |29| ) does not follow the predicted dependence. 
However, for J w 10, one can observe distinct regions of 
(f> with high density of level crossing (these regions are 
distinctly seen in j2(J for J > 6). Upon an increase of J 
these regions narrow down giving the points separating 
the 6-, 8-, and 12-fold configurations. Further increase of 
J leads to a "bunching" of low energetic levels into the 
predicted groups (multiplets) of six, eight, or twelve. 

Not only the numbers of the levels in the multiplets, 
but also the ratios of the spacings between the levels in- 
side the multiplets, the oscillations of the spectra in the 
regime of the multiple tunneling path, and the tunneling 
amplitude in the regime of a single tunneling path obey 
the predictions of our theory. 

For a demonstration we have chosen a set of close val- 
ued J's: J = 23, 47/2, and 24. Figs. Ill 



JT|, and |13| 

for these 



are graphs of the spectra of Hamiltonian 
values of J. The vertical dashed lines are the classical 
boundaries between the different configuration (see the 
diagram of Hamiltonian (||) Fig. ||). A small deviation of 
the dashed line separating the 6- and 8- fold configura- 
tions (06-8 = arctan(3)) towards the 6- fold one is due to 
the fact that, at <f> — 06-8, the depth of the CEF poten- 
tial in the minimum locations of the 6-fold configuration 
is equal to that of the 8-fold configuration. However, the 
intersection of the levels occurs when the ground-state 
energies coincide. See Appendix ^ for details on this 
subject. 

From the pictures one can clearly see the "bunching" 
of the highest and lowest energy levels into the predicted 
multiplets of 6, 8, and 12. The excited multiplets have 
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the same structure which fails only in the vicinity of the 
boundaries between the configurations. The structure 
of the spectra given on Figs. O, O, and 13 looks quite 



similar at this level of "magnification" . To see the subtle 
details predicted in previous sections we should "zoom 
in" the pictures "focusing" on the ground multiplet. 



- 


1 i 
i 

- 1 / ~ 

\ s 1 / 
\/ i >;;'" . 












12-fold 


*^s\ 1 / 

"6-fold 8-pfold 



-2 







Reduced parameter, 
FIG. 11. The spectrum of Hamiltonian (Egh; J ; 



23. 




Reduced parameter, 
FIG. 12. The spectrum of Hamiltonian 



J = 47/2. 
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Reduced parameter, 
FIG. 13. The spectrum of Hamiltonian (pl|; J = 24. 

Firstly, we shift the "center of mass" of the ground mul- 
tiplet to zero (we are not interested in finding the single 
well localization energy) . Secondly, we rescale the shifted 
levels, so that a "visual" comparison of the spacings be- 
tween the levels can be done at different values of the 
reduced parameter (f>. The rescaling is necessary due to 
a large variation of w cx exp(— Jc(<f>)). The calculations 
of the tunneling amplitude for C(0,4) configuration of 
Hamiltonian [| (see Fig. |^) predict a variation of w of 
order 10 6 for J ss 24. The results of this program are 
shown on Figs, [l], and |l| for J = 23, 47/2, and 24 
respectively. The vertical dashed lines (the quasiclassi- 
cal boundaries, see Fig. ||) separate not only the regions 
of different configuration numbers but also the regions 
of the single and multiple tunneling path regimes. The 
regions are enumerated by Roman numerals: I - C(0, 2) 
(single tunneling path regime only), II - C(0,4) single, 
III - C(0,4) multiple, IV - C(0,3) multiple, and V - 
C(0,3) single regimes respectively. Part (a) of each pic- 
ture represents the plot of the rescaling factor which is 
proportional to c(</>); part (b) is the rescaled spectrum. 




-2 

Reduced parameter, 
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FIG. 14. (a) Graph of ln(R)/J, where R is the rescaling 
factor applied to the ground multiplet; (b) Rescaled ground 
state multiplet of Hamiltonian ( p9| ) (the legend shows the de- 
generacies of the levels); J = 23. 
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Reduced parameter, 
FIG. 15. (a) Graph of ln(R)/J, where R is the rescaling 
factor applied to the ground multiplet; (b) Rescaled ground 
state multiplet of Hamiltonian (jj^) (the legend shows the de- 
generacies of the levels); J = 47/2. 




result found in Section || with the numerical calculations 
for J = 24 and 48. The plot is - \n((E 1 - E Q )/A)/J vs. 
u = tan((f>), where Eq and E\ are the energies of the 
ground and first excited states respectively. The differ- 
ence Ei — Eq is 4w according to predictions of Section 
A small discrepancy is due to the coefficient of the expo- 
nential f(u) (w = f(u) exp(— Jc(m))), whose contribution 
decreases oc 1/J. From these data, we can estimate that 
the values of the coefficient f(u) are in a range 0.1 — 3.0. 




Reduced parameter, 
FIG. 16. (a) Graph of ln(R)/J, where R is the rescaling 
factor applied to the ground multiplet; (b) Rescaled ground 
state multiplet of Hamiltonian (^) (the legend shows the de- 
generacies of the levels); J — 24. 



All predictions of Sections |l , [V] and VI (the orderings 
of the levels, the ratios of the level spacings, the oscilla- 
tions of the spectra in some regions, the numbers of the 
oscillations, and the dependence of the scaling parame- 
ter R) find confirmation here. The oscillations are not of 
the periodic form due to a non-trivial (but monotonic) 
dependencies a = a(4>) and £1 = fi(</>). 

A more precise value of c(0) can be easily obtained 
from the single tunneling path part of the spectrum of the 
6-fold configuration. Fig. 17 compares the quasi-classical 



-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 
U 

FIG. 17. Comparison of the quasi-classical and numerical 
tunneling amplitude exponent, u = tan(0). 

All these facts strongly emphasize the validity of the 
developed quasi-classical description of the large spins 
from the theoretical point of view. Now questions arise: 
What is a possible experimental realization? What are 
the limitations of the theory when applied to the real 
systems? We will elaborate these question in the next 
section. 



VIII. EXPERIMENTAL REALIZATION 

A. Feasible experimental systems 

The experimental observation of the predicted effects 
can be done on any system with large values of the angu- 
lar momentum such as rare-earth ions, magnetic clusters, 
or nuclei. The main question is whether the value of J 
is large enough for a given configuration of the external 
field. 

For configuration with small number of minima (2- and 
4- fold configuration), J m 8 satisfies the quasi-classical 
requirement. Such values of J are available, e.g., in rare- 
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earth ions: Dy +3 , Ho +3 , or Er +3 . An example of com- 
pounds with the tetragonal symmetry, where the 4-fold 
configuration is realized, is RENI2B2C, RE stands for a 
rare-earth magnetic element. To suppress the influence 
of the interaction between the magnetic moments the 
magnetic ions should be diluted with similar but non- 
magnetic ones such as La +3 , Lu +3 , or Y+ 3 . The CEF 
effects for this family of compounds were studied in sin- 
gle crystals of Lui_ 2 ,Ho :c Ni2B2C by Cho et al. p3fl . The 
calculated CEF level scheme given there shows that the 
ground state quadruplet is well separated from other ex- 
cited states and corresponds to the multiplet of C(D<4, 2) 
configuration with J = (mod 2) and w » 2K. 

Another family of rare-earth compounds, RESb, of- 
fers the cubic environment. However, it is questionable 
whether the quasi-classical requirement is satisfied since 
even for the highest values of the angular moment ( J = 8 
for Ho +3 ) the multiplicity 2J+1 = 17 is not so large com- 
paratively to the lowest dimension of the cubic configu- 
rations N = 6. The numerical calculations performed in 
the previous section indicate that only for J > 12 there 
are regions of parameter <fi where the 6- and 8-fold con- 
figurations are well defined. To obtain the 12-fold config- 
uration, in the framework of Hamiltonian (|29|), the value 
of J should be increased to about 24. 

Magnetic clusters and molecules offer systems with 
very large total spins and a variety of symmetries. Theo- 
retical calculations (24|] indicate that clusters of 13 atoms 
of transition metals such as Fe, Pd, and Rh may have cu- 
bic symmetry and total magnetic moment of the order of 
Hb per atom. Gadolinium clusters Gd„ (n = 11 — 92) |25| ] 
exhibit large magnetic moments of 0.5-3.0^ per atom 
(which is below the bulk value of 7.55^s but, still offers a 
large value of the total cluster spin) with behaviors rang- 
ing from tight locking to the lattice by crystal anisotropy 
to superparamagnetism (almost free moment). 

Large spins were also observed in artificially grown 
magnetic dots used for observation of the magnetic tun- 
neling 0]. So far, these systems belonged to the lowest 
symmetry class, it is rather tempting to create environ- 
ment of higher symmetry and to use smaller magnetic 
dots like the ones used by I. Schuller and et al ]3G] to 
observe the effects predicted by our theory. 



B. Practicable experiments. Magnetic 
measurements. 

The experimental consequences of the difference 
among the configurations and the spin values can be ob- 
served with many experiments. To name a few these are 
measurements of the spin magnetic moment and mag- 
netic susceptibility, relaxation of the magnetization, elec- 
tron paramagnetic resonance (EPR), and nuclear mag- 
netic resonance (NMR). 

First we discuss measurements of the magnetic suscep- 
tibility. The magnetic susceptibility follows the Curie law 



for temperatures higher than the characteristic splitting 
of the ground-state multiplet Tk B > w (see Appendix [(]). 
Thus, it is (g /ie J) 2 / (k^T) for 1-dimensional configu- 
rations, i.e., C(D N ,N), N = 2,4,6, (gfi B J) 2 /(2fc B T) 
for 2-dimensional ones, i.e., C(Dn,2), N — 4, 6, and 
(g/iB J) 2 /{ik-oT) for the 3-dimensional ones, i.e., for the 
rest of configurations considered in this work. For tem- 
peratures lower than the characteristic splitting T/cb < 
w, the Curie dependence is no longer universal. The non- 
magnetic classes, that is those without magnetic moment 
in the ground state, have their magnetic susceptibility 
saturated to some constant at T — ► 0, whereas, the mag- 
netic ones (with a non-zero magnetization in the ground 
state) still obey the Curie-like behavior. Both the satura- 
tion values and the slopes of the Curie-like curves depend 
upon the configuration of the symmetry group as well 
as upon the equivalence class of the spin. For example, 
the 2-fold configurations C(Djy, N) are non-magnetic for 
J = N/2, N = 2, 4, 6; the saturation values of the mag- 
netic susceptibility are x(T) = (g^B J) 2 /\E(Q)\, where 
E(0) are the corresponding eigenvalues for zero magnetic 
field (see Eqs. (p3)). The magnetic classes of these con- 
figurations, i.e., the ones with J ^ N/2, have the same 
Curie-like dependence: x(T) — (gfJ-BJ) 2 /(ksT). 

In the case of the C(Dn,2) configurations, N = 4, 6, 
the integer spin classes are non-magnetic and the half- 
integer ones are magnetic. The low temperature mag- 
netic susceptibilities of these co nfigu rations can be found 
in the last column of Tables VI, VII for the CiTi^, 2) and 
C(D 6 ,2) configurations respectively. A detailed temper- 
ature dependence of the magnetic susceptibility is shown 
on Fig. n§ for the CiT)^, 2) configuration. 
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FIG. 18. Magnetic susceptibility versus inverse tempera- 
ture. C(Da, 2) configuration. The dotted line is the high 
temperature asymptote. 

The division into the classes of equivalence is more 
subtle? for the high-or der symmetry groups. Ta- 
bles VIII, EX, XI, and XIII collect the low tempera- 
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ture susceptibilities for the C(0, 4), C(0, 3), C(Y , 5) , and 
C(Y, 3) configurations respectively. Figs. |lS| and |20j show 
the details of the transition from the Curie high temper- 
ature regime to the low temperature one for the C(0, 4) 
and C(0, 3) configurations respectively. 
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FIG. 19. Magnetic susceptibility versus inverse tempera- 
ture. C(0,4) configuration. 



0.6 



2.5 



0.5 



pq 
=3. 



0.4 



0.3 



X 0.2 



0.1 







J=3n 
J=3n+l/2 
J=3n+1 
J=3n+3/2 







0.5 



1 1.5 2 

w/(k B T) 

FIG. 20. Magnetic susceptibility versus inverse tempera- 
ture. C(0,3) configuration. 

For magnetic measurements it is important that the 
system is in thermal equilibrium and the range of tem- 
perature T/cb < w is accessible. This requirement means 
that w ~ £oe~ cJ is not too small. On the other hand, J 
must be not less than ~ N/2 to guarantee the validity 
of quasi-classical approximation. At a fixed lower limit 
for experimentally accessible temperature 7} inequality 
N < J < \u(eq / (k-eT)) / c must be satisfied. For rare- 
earth ions £o is the atomic scale of energy and N = 6, 8. 
It gives Ti < £oe~ 4 ~ 100K which is easily satisfied. For 
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Lai-zHo-rN^I^C the estimated numerically w is about 
2K 1 23 1, preliminary experimental results by D. Naugle 
and coworkers give w w IK [ p7fl . 

Gd +3 ion has zero orbital momentum, its anisotropy 
is caused by the relativistic spin-other-orbit interaction 
and corresponding Sq is about 10 time less than the 
atomic scale (eo ~ 1 -f- 10K). The total spin of Gd +3 ion 
S = 7/2, not too large, but may be enough in the case 
of the tetragonal symmetry. The estimated value Ti is 
between 0.1 and IK. 

The anisotropy of a ferromagnetic cluster is induced 
mainly by its boundaries. The anisotropy energy has the 
same magnitude ~ 1 -j- 10K per a site near the boundary. 
For the cluster as a whole this value must be multiplied 
by the number of atoms on one of the faces of the cluster 
(T) which depends on the cluster geometry. An estimate 
can be attained for the series of magic atom- number clus- 
ters @: Ma/-, N = 13, 55, 147, . . .. These clusters are 
obtained by surrounding a core atom progressively with 

additional shells of atoms: Sk = 10fc 2 + 2, k = 1, 2, 3, 

This procedure can be done for icosahedral, decahedral, 
and cuboctahedral packings, which have 20, 15, and, 12 
faces respectively. For M = 55 we find T w 42/16 and 
e « 3 ~ 26K. On the other hand J oc (Af. In the Gd 
cluster £ ~ 0.5 and for N = 55 we find J ~ 27. It is suf- 
ficiently large. The value of w ~ Soe~ cJ with c « 0.3 is 
between 0.001 and 0.01K. For JV = 13, w ranges between 
0.1 and IK. 



C. Spectral analysis. 

The most straightforward experimental approach is the 
spectral analysis. The main difficulty on this way is that 
the scale of the splitting is very different for different 
systems and values of J. Nevertheless, we can expect that 
the spectral frequencies are either in the sub-millimeter 
or in the UHF range. Apart from the direct attenuation 
measurements, it is possible to apply EPR technique. It 
measures the splitting in magnetic field, i.e., magnetic 
moment in some state. The advantage of this method is 
that it does not require too low temperatures. Certainly, 
its sensitivity drops with the growth of temperature, but 
not too fast. 



D. Oscillations of magnetization. 

Let us consider many identical large spins placed into 
external magnetic field along one of the easy directions 
(fe), sufficiently large to polarize them almost to satura- 
tion. If the field is switched off abruptly, each spin re- 
mains in the same state |fc). Since \k) is not a stationary 
state, it will vary in time according to the Schrodinger 
picture: 



\k,t) = J2\ja){ja\k)t 



-iEjt/h 



(30) 
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Here j labels sublevels of one iV-plet and a labels states 
of the jth sublevel. It leads to oscillation of the magnetic 
moment along the k direction in time: 



M(t) = gfx B J ^2 cos 7fcfc' 



(31) 



x ]T (ja\k){j'a'\kr{ja\k')*{fa'\k') 

jaj'a' 



where jkk' is the angle between the directions of clas- 
sical angular momentum in the extrema k and k', and 
ujjji = (Ej — Eji)/Ti is the transition frequency. All spins 
had the same initial state |fc) at the moment when the 
field was switched off, therefore, their magnetic moment 
will rotate coherently creating the macroscopic rotating 
magnetization. Obviously, the rotation energy will dissi- 
pate. Let us estimate the attenuation time r. We assume 
that the spins are embedded into an insulator. Then only 
phonons leads to dissipation. The spin-phonon interac- 
tion energy can be written as follows: 



n 



s _ph — ^u a p J a Jp, 



(32) 



where u a p is the deformation tensor. The value of the 
coupling constant can be estimated as A ~ A/ J 2 , where 
A is the energy difference of two oscillatory levels local- 
ized near one minimum of the potential /(J). A routine 
calculation leads to an estimate of the oscillations life- 
time r: 



hps 5 



(33) 



where p is the mass density of the matrix, s is the sound 
velocity. For typical values p — lOg cm~ 3 , A = 10K, uj — 
w/h = 10 10 s _1 , and s — 10 5 cm s _1 , we find r ~ 10 _1 s. 
The magnetic field must be switched off for a shorter time 
interval. It seems feasible. For Gd we estimated both A 
and w by a factor of 10 smaller than the values we used 
for the above estimate. It gives the attenuation time r 
in the range of few hours. 

In our estimate we assumed that the temperature T is 
less or of the order of w. If it is much larger, the value of 
t ( |33"1 ) must be multiplied by a small factor fou} / {k&T) . At 
a temperature IK with Tilu ~ 0.01 it changes r from few 
hours to a minute, but still leaves this time long. Thus, 
the requirements for temperature is not too restrictive. 

Nevertheless, the observation of the macroscopic oscil- 
lations of magnetization may be obstructed because of 
inhomogeneous line broadening caused by the dipolar in- 
teraction | p8[ . Indeed, the random shift of the frequency 
due to the dipolar interaction is of the order: 



6u> 



g 2 p^J 2 g 2 p 2 i J 2 nx 
hR 3 = h ' 



(34) 



where R is the average distance between large spins, x 
is their concentration per site, n is the density of the 
matrix sites. For g = 2, J = 3.5, and n — 10 22 cm~ 3 we 



find Su) ~ 1.8 ■ 10 10 xs _1 . For x = 0.001 it three orders 
of magnitude less than lo ~ 10 10 s _1 , but it destroys the 
coherence for the time interval 27r(<5w) -1 ~ 10~ 7 s. What 
can be observed after this interval of time is the noise 
in a rather narrow spectral range 5ui given by Eqn. ( [34] ) 
near the frequency to. The noise attenuates during the 
interval r [Eqn.ps. (|3^)] after the pulse of magnetic field. 
Repeating the pulse of magnetic field periodically with 
the period t < r, one can maintain a permanent average 
level of the noise. Also, one can use this narrow-line noise 
to generate a coherent oscillations in a resonator. 



IX. CONCLUSION 

We have shown that large spins (total orbital mo- 
menta) J placed into external fields of high symmetry 
group G display unusual behavior of low-lying and high- 
lying parts of spectra and magnetic susceptibility. These 
parts of spectra are represented by multiplets contain- 
ing N{G,p) states each, where N(G,p) is the doubled 
number of p-fold axes. Each multiplet is splitted into 
sublevels with multiplicities chosen from dimensionali- 
ties of the irreducible representations of the point group 
G and determined by G, J and p. The distances between 
sublevels in the multiplet are proportional to exp (— cJ), 
whereas the distances between multiplets are propor- 
tional to 1/J. The multiplicities at a fixed G and p are 
periodic functions of J with the period p. The relative 
distances between levels are also periodic functions on J, 
but their period is equal to a half of the number of the 
equivalent plaquettes formed by the tunneling trajecto- 
ries and covering the unit sphere. Interesting exclusions 
are the configurations of the octahedron and icosahedron 
groups with p = 2. In this cases the mutual arrange- 
ment of the levels is stochastic, though the multiplicities 
remain fully deterministic. 

In all considered situations with exception of the tetra- 
hedral and hexagonal symmetry with in-plane easy direc- 
tions, the change of large spin J by 1 leads to a drastic 
change in the spectrum and thermodynamic properties. 
We demonstrated that at such a change the magnetic 
susceptibility can cither change its behavior from Curie 
law to saturation or change the coefficient in Curie law. 

Rather special phenomena appear near hyper-surfaces 
in the space of the Hamiltonians which separate regions 
with different configurations of the extrema of the po- 
tential, i.e., regions with different N = N±, N%. Directly 
on these hypersurfaces the number of equivalent extrema 
is equal to N\ + A^. Thus, in a narrow vicinity of the 
hyper-surface there appears a new " class of universality" , 
new set of sublevels with new multiplicities. Moreover, 
we expect a kind of "turbulent" behavior of levels near 
these hyper-surfaces. 

Given the classical Hamiltonian 7Y(J), one can indicate 
a value J c (Tt), starting from which the multiplicities are 
correctly determined by our theory. Though, this value is 
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model-dependent, our numerical calculations show that 
J c ^N(G,p). 

All conclusions of the theory were checked numerically 
for a model Hamiltonian of the cubic symmetry contain- 
ing two invariants (two free parameters) up to J = 60. 
The agreement for the relative distances between the lev- 
els is very good starting from J sw 20. Multiplicities are 
well determined by our theory starting from J k 12 for 
the 6- and 8-fold configurations and from J ss 16 for the 
12-fold configuration. 

We proposed three classes of experimental systems 
which can display the predicted effects. One of them 
is represented by alloys with participation of two lan- 
tanides or actinides, R and R' , so that R has zero orbital 
momentum and its concentration is close to 1, whereas 
the element R' has large J and its concentration is very 
small. In this way the configuration of large spin in a 
symmetric environment is realized. Typical representa- 
tives are Lai-^Ho^N^r^C (tetragonal environment) or 
Lui_ x Dy x Sb (cubic environment). 

The second class of systems are metallic or metallo- 
organic clusters made from ferromagnetic metals. For 
such clusters symmetry can be not only octahedral, but 
also icosahedral, as it is for the cluster Fei3. The clusters 
may have larger total spin than lanthanide and actinidc 
atoms. In both cases we propose to measure spectrum 
of low-lying states (EPR or NMR measurements) and 
also to measure magnetization and magnetic suscepti- 
bility at low temperatures (about 1-2K). Though experi- 
mental difficulties may arise on the way to realization of 
these experiments, we believe that the expected physical 
phenomena are worthwhile to study. 

The third class is magnetic dots used in experiments 
on magnetic tunneling 29 30f| . 

Experimenters should choose optimal values of J to 
ensure the validity of quasiclassical approach: reliable 
separation of the iV-fold multiplets and simultaneously 
not too small values of the tunneling exponent exp(— cJ) 
with c < 0.55 for the cubic symmetry and c < 0.29 for 
the icosahedral symmetry. 

An interesting experimental and may be technical ap- 
plication of our system is the excitation of magnetic oscil- 
lations in a narrow spectral region by pulses of external 
magnetic field. The frequency of these oscillations range 
from 10 7 to 10 n Hz. 
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APPENDIX A: CHARACTERS OF THE MAIN 
REPRESENTATION 

The character of the identity transformation E is triv- 
ial x(E) = dim(W) = N. For a half-integer J, one has 
to find projective representations of the corresponding 
group with the factor set of {±1} or, equivalently, lin- 
ear representations of the group extended by a group 
{E, Q = exp(i2TrJ)E} (the so-called two-valued repre- 
sentations), where Q is a rotation through an angle of 
2tt (here we adopt the notations of Ref . |Ll| ) . Obviously 
X(Q) = cos(2tt J) dim(VF) = qos{2-kJ)N. Other elements 
having non-zero characters are the rotations with respect 
to the axes passing through the directions belonging to 
the set C(G,p). Actually, it is sufficient to consider only 
one element from each class of conjugate elements. To 
calculate the corresponding characters we employ a fol- 
lowing trick. 

Let us consider an element of P which, in our basis, 
corresponds to a rotation with respect to a p-fold axis. 
A set of rotations with respect to this axis forms a cyclic 
subgroup of P, and q is a power of the generator of the 
subgroup C p (minimal non-trivial rotation); q may take 
any integer value. For a given configuration C(G,p) a 
non-zero character may occur only if leaves at least 
two of the extrema i and i unmoved. It means that cither 
the rotation axis passes through i and i, and q ^ pn or the 
rotation is trivial: q = pn, n is an integer. Let us choose 
a tunneling path connecting i and i, and passing through 
intermediate nearest-neighbor extrema i — ► i± — > • • • — > 
in — ► * (see Fig. |2l]). Note, that some of the minima 
may coincide, that is, ij = i% for some pairs j k. The 
rotation transfers each extremum ij (J = 1, 2, . . . , n) 
into i'j leaving i and i unchanged. The two paths form a 
closed loop on the sphere which subtends the solid angle 
of Airq/p. This fact leads to a relation for the oriented 
sum of the phases along the circuit: 



(Al) 



k,ii - <t>i,i' + <t>i n ,% - 4>i'% = J 
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FIG. 21. Transformation of 
i — ► i\ — > ■ ■ ■ — ► i n — > % onto i — > i 
solid angle of the filled area is 4nq/p. 



tunneling path 



The same rotation transforms the phases in a fol- 
lowing way: <j> tj , l]+1 h-> 9^ +1 > 3 = 1, 2, . . . ,71 - 1, 
</>Mi ^ , 4>i n ,i i-> To keep the phases 

unchanged (modulo 27r) and the Hamiltonian invari- 
ant one should apply a gauge transformation hi = 
diag{exp(iai), . . . , exp(aAr)}, then the mappings turn 
into equations: 



3' 3+1 J J ^ 



(A2) 



Substituting these equations into ( |AlD gives: 



ai — oti = J- 



which leads to: 



a io = J h d, a lo 

P 



P 



(A3) 



Mult iply ing hi with exp(— i8) removes the common factor 
8 in ([A5j). Thus, the transformation W = exp(—iS)UC^ 
of the configuration space has only two diagonal elements 
{W}u = {W}^ = exp(i2irJq/p). The character of the 
element W is 2 cos( J2irq/p). All other characters are ze- 
ros since the corresponding rotations do not leave any 
state \k) of the configuration invariant. The characters 
of the representation do not depend upon the gauge: 
Tr(W') = Tr(WWW t ) = Tr(VV). Thus, one can study the 
reduction of the representation W in any specific gauge 
without loss of generality. 



is not known itself. However, the symmetry of the sys- 
tem is of great help and is used throughout the paper to 
assess the connectivity structure. An example, in which 
the geometrically next-nearest-neighbors must be incor- 
porated into the connectivity matrix to geth er with the 
nearest- neighbors is discussed in Section V B . 

The connectivity matrix of C(0,4) configuration is 
quite simple since only the geometric nearest-neighbors 
should be connected. For the minima enumeration given 
on Fig. |] the connectivity matrix is: 



Oil 

i 6 i 
i i 6 








i i 





, i = 


i i 



(Bl) 



This is, actually, the Hamiltonian for J = (mod 4), 
without a multiplier of vo, since the phase structure is 
absent or of no importance for the respective cases. For 
all other J's one should find the twelve phases of the tun- 
neling amplitudes. Five of the phases may be set to zero 
due to the gauge freedom with the only constraint that 
Eq. ( |l2| ) must be satisfied. In our sample case, we null the 
following ones: 0i,3, 4>i,i, 4>i,5, <fii,6, 4>2,5- The rest of the 
phases is obtained from the seven independent plaquettes 
(each plaquette gives an equation of type of Eq. (p"2|)); 
these are: 02,3 = — ttJ, ^ 2 ,4 = ttJ, </>2,6 = 2ttJ, </> 3i5 = 
ttJ/2, <j>3 t e = — 7rJ/2, 04,5 = — 7rJ/2, 04 ; 6 = ttJ/2. Thus, 
the defined Hamiltonian is: 



H? 









a Z7T J 

—inj 



1 e" l27rJ 



1 





a iTrJ/2 



1 




e iTr.//2 
-ittJ/2 



1 
1 

e in,J/2 
e -tirJ/2 e i7rJ/2 



1 

i27T J 

-inJ/2 



(B2) 



APPENDIX B: AN EXAMPLE OF A 
HAMILTONIAN FOR C(0,4) CONFIGURATION 
AND CALCULATION OF ITS SPECTRA. 

In this appendix we consider a detailed construction of 
the reduced Hamiltonian for C(0,4) configuration. First 
we define the connectivity matrix of the system deter- 
mining decide which minima ought to be connected by 
tunneling paths. This is usually done by connecting the 
nearest-neighbor minima. In some cases, however, the 
geometric closeness on the sphere is not a good crite- 
rion for connecting minima. A "sure-fire" criterion is the 
path integral approach which determines the amplitude 
of a spin transition from one localized state to another 
by summation of the contributions of all trajectories con- 
necting the minima. This technique gives the exact solu- 
tion of the problem, but it is very complicated. We could 
use instead a semiclassical method of finding trajectories 
with minimal imaginary classical action, but the action 



Finally, one needs to find the eigenvalues of the Hamil- 
tonian. The diagonalization can be performed by any 
symbolic solving system or even "manually" since the 
Hamiltonian can be factorized. Also, a trick of a purely 
geometric origin can be used: if one views Hamiltonian 
H (w = 1 and \(H)ij\ = 1 if (W)y ^ 0) as a weighted 
connectivity matrix of a graph, then: 



]T 9i E? = TrH n = 2J2J2 cos (°^ ) = J " 



(B3) 



where n = 0, 1, 2, . . ., Ei is the ith distinct eigenvalue 
of TL of multiplicity gi, the first sum on right-hand side 
is over the vertices of the graph, the second is over all 
closed loops of n walks running through vertex j and 
is the flux passing through the fcjth loop. If (H)u ^ 0, 
i = 1, . . . , N, some extra weights nee d to be applied for 
each loop. An advantage of formula (|B3| ) is that /„ are 
gauge invariant s ince all Q&. are gauge invariant. Applied 
to Hamiltonian (BS) for n < 4, the formula yields: 
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i i 

J2^ E i = 24,^.9^f = 48cos(7rJ/2). 



Let us find the spectru m of H amiltonian (B2) for J 



0. We know from Section 

1, g 2 = 2, and g 3 = 3. 



[VB1 



Too, 



of Section IV that g\ = 
it is clear that E\ = 4 



(sum of the matrix elements in each row or column of 
Hamiltonian (B2) is equal to 4 for J = 0). Two equations 
for E2 and £3 are: 



4 + 2E 2 + 3£ 3 = 0, 16 + 2££ 



3E£ = 24. (B5) 



Solving (|B5| ) and checking the roots against the last equa- 
tion of (|B4| ) gives _E 2 = —2 and E 3 = 0. For J = 2, the 
spectrum is inverted (see Section |v|): i?i = —4, E 2 = 0, 
and £3 = 2. For all other J's the solution is trivial since 
there are only two distinct eigenvalues. 



APPENDIX C: LOW TEMPERATURE 
MAGNETIC SUSCEPTIBILITY. 

For high temperatures k B T > w, the moments are 
purely classical and show Curie magnetic susceptibility: 



Xc = {p.BgJ) 2 /(dk B T). 



(CI) 



where d is the dimensionality of the system, fi B is the 
Bohr magneton, k B is the Boltzmann constant, and g is 
the gyro-magnetic ratio. For low temperatures k B T < w, 
the quantum effects change the response drastically. The 
susceptibility saturates for the classes without magnetic 
moment in the ground state to the value: 



Xs = 



90 



(C2) 



where go is the degeneracy of the ground state and Xi is 
the susceptibility of the ith member of the ground-state 
multiplet. For the classes with the magnetic moment in 
the ground state Curie susceptibility persists at k B T < 
w, but its slope is different: 



Xi 



-, 90 

-E™. 2 /(feT), 



(C3) 



where rrij is the moment of the ith member of the ground- 
state multiplet. 



-J/04 



H? 3 



-a — b 
9 



- (2a + 3b)(x 2 
-{Ab-a)(x 2 



(Dl) 
(D2) 



where x and y are local Cartesian coordinates. Treat- 
ing these terms as the effective potential energy of the 
quantum-mechanical problem, one can identify it with 
that for a 2-dimensional isotropic harmonic oscillator 
(the kinetic energy is due to the Wess-Zumino term |2(| 
or Berry phase, its exact form is of no importance here; it 



suffices to know that this term is identical for both (Dl) 
and (D2), thus, providing identical effective masses M). 
The potential energies of the two harmonic oscillator 
problems are equal on a line b = 3a. The squares of the 
effective frequencies are 22a/M and 88a/(3M) (a > 0), 
respectively for C(0, 4) and C(0, 3), on this line. Hence, 
the energy of the ground as well as the spacing between 
successive levels is larger for C(0, 3) configuration on this 
line. These arguments qualitatively explain the deviation 
of the line b = 3a from the point of the level crossing to- 
wards the 6-fold coordination region. 

The boundary of the transition from one configuration 
to another, (a line in our 2-dimensional a-b space) marks 
a singularity in the "flow" of the level multiplicities across 
the parameter space. The levels of two different coordi- 
nations must match exactly at this surface. We find the 
spectra of this intermediate configuration in the case of 
N = N(0, 4) + N(0, 3) = 14. 

The surface of the sphere is covered with twelve con- 
gruent even-sided plaquettes. Hence, the period of spec- 
tra is J = 6 and all spectra are symmetric (see Sec- 
tion |v|). The spe ctra of the 14- fold configuration are 
collected in Table. XVI for the non-equivalent J's. 



TABLE XVI. The spectra of a "hybrid" C(0, 4) + C(0, 3) 
configuration. 

J Eigenvalues(degeneracies) 

±2^™ (1, 1), ±2w (3, 3), (6) 

1 ±(l + v%)w(3,3), ±(1- y/S)w(3,3), 0(2) 

2 ±V6w (2, 2), ±2w (3, 3), (4) 

3 ±2w ( 6,6), 0(2 ) 

1/2 ±wy/6 + 2^(2, 2), ±w V / 3 _ V5 (4, 4), (2) 
3/2 ±wy / 6(4,4), (6) 

5/2 ±1^6-2^(2, 2), ±w\/3 + v / 3(4, 4), (2) 



APPENDIX D: ON THE INTERSECTION OF 
THE MULTIPLETS OF C(0,4) AND C(0,3). 

In a close vicinity of a minimum of C(0, 4) and C(0, 3) 
configurations, Hamiltonian (|^) has the following form: 
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